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Low temperature properties of the spin- 1/2 frustrated ladder with ferromagnetic rungs and legs, and 
two different antiferromagnetic next nearest neighbor interaction are investigated using the modified spin 
wave approximation in the region with ferromagnetic ground state. The temperature dependence of the 
magnetic susceptibility and magnetic structure factors is calculated. The results are consistent with the 
numerical exact diagonalization results in the intermediate temperature range. Below this temperature 
range, the finite size effect is significant in the numerical diagonalization results, while the modified spin 
wave approximation gives more reliable results. The low temperature properties near the limit of the stability 
of the ferromagnetic ground state are also discussed. 
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1. Introduction 

Among a variety of models and materials with strong 
quantum fluctuation, quantum spin ladders have been 
attracting the interest of many condensed matter physi- 
cists. Recently, a spin-1/2 ferromagnetic ladder com- 
pound Cu n Cl(0-mi)2(/i-Cl)2 (mi — 2-methylisothiazol- 
3(2iJ)-one) is synthesized with possible frustrated inter- 
rung interaction. 1 

Although frustrated and unfrustrated spin ladders 
with nonmagnetic ground states have been extensively 
studied, 2-5 those with ferromagnetic ground states have 
been less studied. This would be due to the simplicity of 
the ferromagnetic ground state. Even if the ground state 
is ferromagnetic, however, the excited states and finite 
temperature properties should be influenced by frustra- 
tion. Especially, near the stability limit of the ferromag- 
netic ground state, we expect characteristic temperature 
dependence of physical quantities. In the present work, 
we investigate the finite temperature properties of frus- 
trated spin ladders with ferromagnetic ground states us- 
ing the modified spin wave (MSW) approximation 6-11 
and the numerical exact diagonalization (ED) calcula- 
tion for short ladders. 

The MSW approximation was first proposed by Taka- 
hashi 6 to investigate the low temperature properties of 
ferromagnetic chains. If the conventional spin wave ap- 
proximation is applied to one-dimensional ferromagncts, 
the thermal fluctuation of magnetization diverges at fi- 
nite temperatures. This is natural, considering the ab- 
sence of finite temperature long range order in one di- 
mension. 12 This divergence, however, prevents the cal- 
culation of physical properties such as magnetic suscep- 
tibility and magnetic structure factor at finite tempera- 
tures. To circumvent this difficulty, Takahashi proposed 
to introduce the chemical potential /i for magnons and 
to fix /i by imposing the constraint that the expectation 



Fig. 1. Structure of frustrated ladder studied in the present 
work. 



value of the total magnetization S* ot vanishes as 







(1) 



taking account of the Mermin- Wagner theorem. 12 
This procedure is quite successful for one and two- 
dimensional ferromagnets 6-8 and two-dimensional anti- 
ferromangets. 9-11 Recently, this method has been suc- 
cessfully applied to dimerized ferromagnetic chains. 13 In 
the present work, we employ this method to calculate the 
low temperature magnetic susceptibility of the present 
model with ferromagnetic ground states. We also calcu- 
late the magnetic structure factor to confirm that the 
short range correlation is properly described within the 
MSW approximation. 

This paper is organized as follows. In the next section, 
the model Hamiltonian is introduced. The MSW analysis 
is explained in §3. In §4, the MSW equations are numer- 
ically solved and the results are compared with the ED 
calculation for short chains. The last section is devoted 
to summary and discussion. The low temperature expan- 
sion of the MSW equations is explained in Appendix. 
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2. Hamiltonian 

We consider the frustrated ferromagnetic ladder with 
spin- 1/2 described by the Hamiltonian 

L 
i=l 

+ JiSi,lSi,2 + JdlSi,lSi+l,2 + "^d2'S'i+l, 1^1,2] , (2) 

where Sij (i = l,...,L,j = 1,2) are the spin opera- 
tors with magnitude S. For the numerical calculation, 
we only consider the case of S = 1/2. The number of 
the unit cells is denoted by L. In this paper, we focus on 
the case J\, J r < (ferromagnetic) and Jdi, Jd2 > (an- 
tiferromagnetic). The lattice structure is schematically 
depicted in Fig. 1. We also use the parameterization 

Jdi = Jd(l + 5), Jd2 = Ml-S) (3) 
if convenient. 

3. Modified Spin Wave Analysis 

3. 1 Formulation 

We employ the standard Holstein-Primakoff transfor- 
mation, 



S+.=^.+^.=^/ i;j (S)a i , j , 
Sr.=Sf tj -iSf tj = y/2Say id {S), 



kj(S) = V 1 - (^J-^L-Oij (i - h-,L;j = 1,2), 

(4) 

where aj ■ and dij are magnon creation and annihilation 
operators on the i-th rung and the j-th leg. After the 
Fourier transformation with respect to i, the Hamilto- 
nian is rewritten as 



H — Ho + Hi i 



^o= 5(4,4,2) (| t)(: ; ; 



(5) 
(6) 



Hi 



= -izEEK l(fe " 9) (i- e 4fe ')(i- e ") 



/c/c'g j — 1 



+e lk (l - e - 4(fc+9) )(l - e -* (fe '- 9) )} 
x a l+ q ,j a k>- g ,j a k',jak,j 

~ E {4+g,l4'- 9 ,2 

X(flfc,l - flfc,2)(flfc',l - flfc',2) + h.c} 



-*(*'-«)«t -t 



a fc+ 9 ,l a fc'-9,2 



kk'q 



x (a fc ,i - e lfc a fcj2 )(afcM - e lfc 'a fe ',2) + h.cj 



-i(fc'-<,) t a t 



x (a fc , 2 - e* fc a M )(a^2 - e ife 'a fe / ; i) + h.cj 



(7) 



up to the quartic order in a k j and a\ i . Here, i?i and E 2 
are defined by 



E x = {-2J\(l - cos fc) - J r - J d i - J d2 } 5, 
£ 2 =(J r + e ifc J dl + e" ife J d2 )S. 
We introduce the unitary transformation 

vt 



(8) 



1 



a l 



aLe 4 2 + 



t 

a[ 2 e ^ J , 

)■ 



(9) 



(10) 



and the corresponding magnon number operators 

h ak = a\a k , hp k = p\f3 k . 

The phase </>(fc) is determined afterwards to minimize 
the free energy. We assume the density matrix for the 
magnons in the form of a product of single-magnon den- 
sity matrices as 

P= ^ H- P afe("afe)-P/3fc(«/3fc)|{«afc,«/3fc})({^afe,«,3fc}| 
{n ak ,n pk } k 

(11) 

where 

|{n afe ,n^}) = nKfe!^!)-'(4) nafc (^) n ' 3 ' t |0) (12) 
k 

is the eigenstate of n ak and hp k with eigenvalues n ak 
and np k . In the following, the expectation value with 
respect to (11) is denoted by (...). The probability that 
the state with wave number fc is occupied by n bosons is 
denoted by P vk (n) (y = a or /3). Therefore, the following 
normalization conditions are required for all fc. 



J2P ak (n) = l , ^P /8fc (n) = l. 

n=0 n=0 

Accordingly, the free energy is given by 



(13) 



F=(Ho) + (H 1 )+T ^ J2J2P,k(n)^P vk (n) 

v=a,j3 k n=0 

(14) 

where the expectation values ("Ho) and (Hi) are ex- 
pressed as follows 

(Ho) = + \{E-2e- l ^ k) + E* 2 e^)5h- k } , 

(15) 

Y (h k - <5n Qfe cos0(fc)) 

k 

^ (n k - 5n k cos (<j>(k) - fc)) 



4L 

Jdl 

1l 



k+q,2 u k'-q,l 



kk'a 
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(16) Substituting (21) and (22) into (24) and (25), we find 



n vk = [n. 



^ (n k - Sn k cos (<f>(k) + k)) 

<Sintra(<7) = S + ^ ^2 n k+q/2 h k - q / 2 , (26) 

" fc> = ? ^ fc(n) ( " = ^ (17) WW - + ^fn ter ( 5 ), (27) 

We also define gc^^ = J_ £ cos + g/2) _ _ g/2)) 

«fc = n ak + hp k , Sh k = h ak - np k (18) k 
for convenience. The condition (1) reduces to x Sn k+q/20~ n k- g /2, (28) 



S 2L 

k 



l -^k. (19) S? zltel (q) = ^ E Y, S m(<f>(k + q/2)-<P(k-q/2)) 

' U k 



Introducing the Lagrangian multipliers fi^k an d M which x Sn k+q / 2 Shk- q /2- (29) 

account for the constraint (13) and (19), we minimize the ft should bc nQtcd ^ hag an imaginary part 

following quantity W with respect to P vk and 0(fc). becauge of ^ &bsence rf space inversion symmetry x ^ 



—a;. 



w = p- ^ 5>^E p " fe ( n )-<" E EE^w- 

v=a£ k n=o v=a,p k n=o 3^ Lowest order approximation 

(20) 

To the lowest order approximation, we only consider 

the Hamiltonian Ho and impose the constraint (1). In 
3.2 Spm-Spm Correlation Function and Magnetic Sus- ^ foUowing) we call this ap p rox i mat i on the M SW0 ap- 

ceptibility ■ ,. 

* y proximation. Minimizing 

The spin-spin correlation function is expressed in 

oo 

terms of the magnon occupation numbers as TT7 , V* n c \ 

(SnSi'i) =(Si2Si>2) = ^^-^cosfc(x;< -x 4 )n fc ^ , 

(21) 



v=a,B k n=0 



^ E EE nP ^( n ) ( 3 °) 

y=a,/3 fe n=0 



(SuSm) = ( ttf E cos (^( fc ) ~ fc ( 
V k 



x v - Xi)) Sn k 



with respect to <f), we find 
dW d (Ho) i 



(£ 2 e-^( fe ) - £ 2 V* (fc) )c5n fc = 0. 



(22) 50(A) 50(A) 2* 

(31) 

where Xi is the position of the i-th site. Using these ex- 
pressions, the magnetic susceptibility per spin \ IS ex_ This ^ ea( ^ s to 

pressed as Jr gin ^ + sin ^( fc ) _ fc ) + j d2 sin (0( fc ) + fc) = o, 

x = (.9Mb) 2 y y y ^ g , ^ (32) 

2PP 4^ / — f iJ * J which determines the phase 0(fc) as 

»»' j=i j — l ~ \ / 

j^x± ±<!„ . s „„, «-«*> - SXVZL ^- < * < *«■ 

n> j=ij'=i (33) 

(ou B ) 2 99 , ^ Minimizing W with respect to P^fc, we find 
- ■ 2^("afc + n^fc + + "/*)> (23) 



ne v (k) + T(lnP uk + 1) - n/i - = (f = a,/3), 



(34) 



where is the Bohr magneton and g is the gyromag- 
netic ratio. 

where 

3.3 Magnetic Structure Factor 1 d (Ho) 

We define the intraleg and interleg magnetic structure £ v{k) — — (35) 

factors as 

are two branches of the magnon excitation energy given 
Sintrafa) = T E ( S o,i s i,i) exp(ix^) by 



e» (fc) = Pi ± i(P 2 e-^( fe ' + P*e^) 



= tE (So.2Si. 2 ) exp(ixiq), (24) 

L V = S{-2Ji(l - cosfc) - J r (l=Fcos0(fc)) 

5 intcr (<z) = I V (S 0l iS iia ) exp^g). (25) " Jdl ^ T COS ( fc " " Jd2 ^ T C0S ( fc + 
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The Hamiltonian Ho is rewritten as 

k 



(a fc /3 fc ). (37) fc ? 



Eliminating [i vk from (34) using the condition (13), we 
find 

P vk {n) = {l - e -(^(fc)-M)/r J e -n( e „(fc)-*0/r (38) 
The expectation value of n„fc is given by 



i v k = ^nP uk (n) 



1 



71=0 



e (e„(*0-*O/T - 1' 



(39) 



The chemical potential is determined using the condi- 
tion (19). 

3.5 0(S°) approximation 

In this approximation, we include the terms of 0(S°) 
in the Hamiltonian i.e. Hi. We call this approximation 
the MSW1 approximation. Minimizing 

oo 

Wi = <Wo> + <^i>-T5- ^ $>"*£ P "*( n ) 

i/=a,/3 fe n— 

oo 

E EE nP ^( n ) ( 4 °) 

i/= a,/3 fe n— 

with respect to </>(&), we find 

awi _ d (Ho) _ 

d<t>{k) d<t>{k) d<t>{k) K ' 

which yields 

J r sin (j)(k)S' 2 + J dl sm(<j>(k) - k)S' 3 

+ J d2 sin (4>(k) + k)S' 4 = 0, (42) 

where 

Si = -^-j- ^2 cos kn k , 



S' 2 = -^-^2cos((j)(k))Sfi k , 

k 

s 'z = E c ° s (^( fc ) - fc )^fc' 
,s 4 = ^-E cos W fc ) + fc )^ fe - ( 43 ) 



2i 

k 

From eq. (42), </>(fc) is given as 

{JaiS'z - J d2 S' 4 )s'mk 



tan <j)(k) 



(j T s 2 + (j dl s 3 + j d2 s' 4 ) cos ky 

Minimizing W\ with respect to P vk , we find 
dW x 



(44) 



where 



ni vk + T(ln + 1) - n/i - = 0, (45) 



< = - *J,, ((Ko) + (Hi)) (46) 



is the dressed single magnon excitation energy given by 
so, (k) = -2Ji(l - cos k)S[ - J r (l T cos(0(fc))S , 2 



- J dl (l T cos(0(fc) - fc))5g - J d2 (l T cos(0(fc) + fc))S£. 

(47) 

Expressions for P vk (n) and n vk are obtained by replacing 
£vk by ffj/fe in (38) and (39), respectively. The phase 4>(k) 
and chemical potential \i are determined by solving (42) 
under the condition (19) numerically. 

3.6 Low Temperature Behavior 

We employ the MSWO approximation to investigate 
the low temperature behavior. Among two branches of 
the magnon excitation, £p(k) has an energy gap at k = 0. 
Therefore, we only consider e a (k) at low temperatures. 
Near k = 0, the dispersion relation is given by 



e a (k) ~fSk 2 



where 



Ji 



Jr(Jdl + Jd2) +4J d iJ d 2 



(48) 



(49) 



2(J r + Jdl + Jd2) 

is the effective ferromagnetic exchange interaction be- 
tween the ferromagnetic rung dimers consisting of S^i 
and Si_ 2 . For negative (ferromagnetic) J\ and J r , ^ is 
positive as long as the frustrating antiferromagnetic in- 
teraction Jdi and Jd2 are small. However, with the in- 
crease of Jdi and Jd2, ^ decreases and vanishes at 

J r (Jdi + Jd 2 ) + 4J d iJd2 J r + 2J d (l-<^ 



J] 



2( J r + Jdl + Jd2) 



J r + 2 Jd 



-Jd- 

(50) 



This point is the limit of the stability of the ferromagnetic 
ground state. 

For ^ > 0, the low temperature expansion can be car- 
ried out in the same manner as the ferromagnetic chain. 
The details are explained in Appendix. As a result, we 
obtain 



8(. 9M B) 2 Sy L 3 C(|) / T | 3 C 2 (|) T 



3T 2 



45 y JS 16S 2 7T JS 

(51) 



Namely, the susceptibility is proportional to T~ 2 at low 
temperatures. Defining J? = ^/2,S e s = 2S, the sus- 
ceptibility per rung is given by 



Xrung 2% 



{g^?Sj s J f 3C(j) 



T 



3C 2 (|) 



3T 2 



727 \/ 2/5 e 3 ff 2^ 2 JS\ S/ 

(52) 



This coincides with the low temperature susceptibility 
of the ferromagnetic chain with spin S e s and effective 
exchange interaction ^ . 6 



n dP uk (n) 
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xT 2 [cm" 3 mol K 2 ] 



J--14.47K J,=-5.64K 
J d1 =3.91K J d2 =0 
g=2.13 




□ MSWO 
■ MSW1 

O Extrapolation from L=6,8,10 
• Experiment 



1 



2 T 1/2 [K 1/2 ] 3 



Fig. 2. Magnetic susceptibility for exchange parameters corre- 
sponding to [Cu n Cl(0-rm)2 0-Cl) 2 ] calculated by MSWO (open 
squares) and MSW1 (filled squares) approximations. Open circles 
arc Shanks extrapolation from the ED data for L = 6, 8 and 10. 
The filled circles arc experimental data taken from ref. 1. 



4. Numerical Results 

4-1 Results for the parameter set corresponding to 
Cu ll Cl(0-mi) 2 (v-Cl) 2 

The magnetic susceptibility is calculated for the pa- 
rameter set determined in ref. 1 for Cu n Cl(0-mi)2(/z- 
Cl)2 as shown in Fig. 2. The susceptibility calculated 
by two MSW approximations and the Shanks extrap- 
olation 14 from the ED data for L = 6, 8 and 10 are 
shown. The experimental data for Cu n Cl(0-7m) 2 (£«-Cl)2 
are also shown. Although Fig. 2 contains only a limited 
number of experimental data, there are many data points 
at higher temperatures and the exchange constants are 
determined by fitting them with the ED calculation as 
explained in ref. 1. 

The MSWO approximation reproduces the overall tem- 
perature dependence including relatively higher temper- 
ature regime. At low temperatures, where the ED results 
are strongly size dependent, both MSW approximations 
give the consistent results. The result of the MSW1 ap- 
proximation is more smoothly connected to the ED data. 
However, it strongly deviates from the ED data as the 
temperature is raised. Actually, the MSW1 equations 



have no solution other than 



0,S[ 



S' 2 



S'=0 for 



T > T* ~ AK. This implies a phase transition at T = T*. 
However, this transition is spurious because the present 
system is one-dimensional and no phase transition should 
take place at finite temperatures. Therefore, the MSW1 
approximation is unreliable at T>0(T*). This spurious 
transition takes place even in the absence of frustration 
within the MSW1 approximation. A similar kind of spuri- 
ous transition takes place in the square lattice Heisenberg 
antifcrromagnct. 9 Hence, this is an artifact of the MSW1 
approximation rather than the effect of frustration. On 
the contrary, the MSWO approximation predicts no phase 
transition and gives the results qualitatively consistent 
with the ED results up to relatively high temperatures. 
To confirm that the magnetic short range correlation 



(a) ■ ED L=6 

□ ED L=8 
O ED L=10 



O MSWO 
• MSW1 

J--14.47K J,=-5.64K 
J d1 =3.91K J d2 =0 
T=2K 




(b) ■ ED L=6 

□ ED L=8 
O ED L=10 



O MSWO 
• MSW1 

J--14.47K J,=-5.64K 
J d1 =3.91K J d2 =0 
T=3K 




(c) ■ ED L=6 

□ ED L=8 
O ED L=10 



O MSWO 

J r =-14.47K J,=-5.64K 
J d1 =3.91K J d2 =0 
T=4K 




Fig. 3. Magnetic structure factors with exchange parameters 
corresponding to Cu n Cl(0-mi)2 (/i-Cl) 2 calculated by MSWO 
(small open circles) and MSW1 (small filled circles) approxima- 
tions. The ED data for L = 6 (big filled squares), 8 (big open 
squares), and 10 (big open circles)are also shown. Temperatures 
arc (a) 2K, (b) 3K, and (c) 4K. 



is appropriately taken into account by the MSW approxi- 
mations, the magnetic structure factors are shown in Fig. 
3. With the decrease of temperature, the short range 
ferromagnetic order develops as expected. The results 
of ED calculation are also presented in Fig. 3 for the 
same choice of exchange constants. Although the MSW1 
approximation reproduces the ED results better than 
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Ji 



-2 



(a) 8=0.8 J d =1 


nonmagnetic 


Ferro 


\ **• 


-6 




(b) 5=0.5 J d =1 


nonmagnetic 


Ferro 


\ \ 


6 


" 4 



Fig. 4. Ground state phase diagram for = 1 with (a) 5 = 0.8 
and (b) <5 = 0.5. The thick line is the ferromagnetic-nonmagnetic 
transition line and the dotted line is the limit of the stability of 
the ferromagnetic state given by (50). 



MSW0 at T = 2K, it becomes poor at T = 3K where the 
susceptibility also deviates from the ED result. On the 
other hand, the MSW0 approximation gives reasonable 
agreement with the ED results even at T = 4K where 
the MSW1 equations only have an unphysical solution. 

4-2 Results at the ground state phase transition point 

With the increase of the frustration, the ferromagnetic 
ground state becomes unstable against magnon creation 
at J\ = J\ s and the transition to the nonmagnetic ground 
state takes place. However, the ground state phase tran- 
sition does not always take place at J\ = J\ s . The ground 
state can change from ferromagnetic to nonmagnetic at 
J\ — Ji c (< Ji s ) where the ground state energy of the 
nonmagnetic state becomes equal to that of the ferro- 
magnetic one. Examples of ground state phase bound- 
aries and stability limits of the ferromagnetic state are 
presented in Fig. 4(a) for 8 = 0.8 and in Fig. 4(b) for 
5 = 0.5. The ferromagnetic- nonmagnetic phase bound- 
ary is determined by the ED method for L — 10. It is 
checked that the size dependence is negligible in the scale 
of this figure. 

First, let us examine the behavior of the magnetic sus- 
ceptibility at J] = Ji c in the case J\ c — J\ s . On the 
line J\ — Ji s , the excitation energy e a (k) is proportional 




Fig. 5. Magnetic susceptibility at the ferromagnetic- 
nonmagnetic transition point (J r = — 8, <5 = 0.8, J\ = J\ C (J T ,5) ~ 
— 1.213333) which coincides with the limit of the stability of the 
ferromagnetic ground state J\ s . The MSW0 and MSW1 results are 
shown by open circles and double circles, respectively. The extrap- 
olated values from the ED results with L = 6, 8 and 10 are shown 
by open squares. The curves fitted by \ ~ (A + BT 1 / 4 )T — 4 / 3 
are also shown. For the fitting of the MSW0 data, A is fixed to 
0.1795937 as obtained in Appendix. 



MSW0 approximation behaves as T -4 / 3 as shown in Ap- 
pendix. As an example of the case J\ c = J\ s , the tempera- 
ture dependence of the susceptibility for J r = —8, Jd = 1 
and 5 = 0.8 and J\ = J\ s ~ —1.213333 is shown in Fig. 5. 
The results calculated by the MSW0 and MSW1 approx- 
imations, and those extrapolated from the ED results for 
L = 6, 8 and 10 are shown. All results show the behavior 
X ~ T~ 4 / 3 . However, the ampltude of the susceptibility 
is overestimated by the MSW approximations. We may 
understand this discrepancy in the following way: 

Even at J\ = J\ s , the ferromagnetic state remains 
one of the ground states. In addition, the nonmagnetic 
state which replace the ferromagnetic states for J\ > J\ s 
comes into play. However, this state has no magnetic 
moment and does not contribute to the susceptibility. 
Also, the low-lying excited states around the nonmag- 
netic states have small magnetic moments and do not 
have significant contribution to the susceptibility. Never- 
theless, these states have finite statistical weight. There- 
fore, the contribution from the ferromagnetic state and 
the excitations around it, which is correctly described 
by the MSW approximations, can reproduce the lead- 
ing temperature dependence of magnetic susceptibility, 
while its actual amplitude is reduced from the results of 
the MSW approximations. 

These nonmagnetic states have antiferromagnetic or 
incommensurate short range correlations induced by 
frustration. To get more insight into their effects, we cal- 
culate the magnetic structure factor for this parameter 
set as shown in Fig. 6. It is clearly observed that the 
q = component, which is responsible for the suscep- 
tibility, is overestimated in the MSW0 approximation. 
It should be noted that the size dependence of the ED 
results for S'intra(<Z = 0) and Si ntC r(q = 0) is almost negli- 
gible. Therefore, this discrepancy is not attributed to the 
finite size effect. On the other hand, the large q compo- 
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-8 J d =1 5=0.8 
Ji=Ji c (8,J r ) T _ 
=-1.213333 



L=10 8 6 

- O O o 
T=0.2 — □ □ □ 




0.5 



+ MSWO 
□ ■ L=10 
O • L=8 

A A L=6 



T 1/3 C MSW0 

T S in , rj ,(q=0) 



intraV 
+ 



fi * ^ 



D t-1/3 ED 



(q=0) 



T 1/3 AS intra (q=0) 



0.2 



0.4 



Fig. 6. Magnetic structure factors on the ground state phase 
boundary calculated by MSWO approximation(lincs) and ED(opcn 
symbols) methods. Temperatures are T = 1 (dotted lines, circles), 
0.2 (broken lines, squares) and 0.1 (solid lines, triangles). The big, 
medium and small symbols represent the system sizes L = 10, 8 
and 6, respectively 



Fig. 7. Temperature dependence of Si n tra(? = 0) on the ground 
state phase boundary calculated by the MSWO approximation (+) 
and ED methods (open symbols). The difference ASi ntra (<? = 0)(= 
S^™(q = 0) - S^Jq = 0)) between the MSWO and ED results 
is also shown (filled symbols). All quantities are multiplied by T 1 / 3 . 



nents are underestimated in the MSWO approximation. 
Considering the sum rule 

SintraG?) = (5 ,l5 ,l) - S(S + 1), (53) 

<i 

the enhancement of q ^ components suppresses the 
q = component for S ln t ia (q). F° r Winter (q), we have 

$>mtcr(<Z) = (S ,lSo,2>. (54) 
Q 

In this case, lhs is not a constant. However, this should 
be close to 1/4 if J r is strongly ferromagnetic as in 
the present case. Hence, the similar explanation is nat- 
urally valid for Si nte r(q)- Thus, we may conclude that 
this discrepancy comes from the frustration induced an- 
tifcrromagnetic or incommensurate short range order 
which is not appropriately described by the MSWO ap- 
proximation. Fig. 7 shows the temperature dependence 
of Sintra(<Z = 0) calculated by the MSWO approxima- 
tion (Sf^°(q = 0)) and that calculated with the ED 
method (S;™ a (<7 = 0)). The difference AS intra (g = 0)(= 
Snfrf °(q = 0) - S™ a (g = )) is also plotted. All 
quantities are multiplied by T 1 / 3 . This plot shows that 
T 1 / 3 A5 l i ntra (<7 — 0) does not diverge in the low temper- 
ature limit. Hence, the correction to Sf^™°(q = 0) does 
not diverge with the power stronger than T -1 / 3 in the 
low temperature limit. This means that the power of the 
leading term of \ <~ T~ 4 / 3 is unaffected by this correc- 
tion. 

If the limit of stability of the ferromagnetic phase J\ s 
is away from the ferromagnetic-nonmagnetic transition 
point J\ c , the susceptibility behaves as T~ 2 as in the case 
of ferromagnetic ground state. An example is shown in 
Fig. 8 for J r = -4, J d = 1 and S = 0.5 and J\ = J lc ~ 
-1.41. 

5. Summary and Discussion 

Low temperature properties of the frustrated ferro- 
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Fig. 8. Magnetic susceptibility at the ferromagnetic- 
nonmagnetic transition point (J r = — 4, J\ = —1.41,5 = 0.5 
) which is away from the limit of the stability of the ferromagnetic 
ground state. The MSWO results are shown by open circles and 
extrapolated values from the ED results with L = 6, 8 and 10 
are shown by open squares. The solid line is the fitted line by 
X oc const. X T~ 2 



the MSW and ED methods. The magnetic susceptibility 
and static magnetic structure factors are calculated. It is 
found that the MSW method gives reasonable agreement 
with the ED calculation in the intermediate temperature 
regime below which the ED results show considerable size 
dependence. On the contrary, the MSW method becomes 
more reliable in the low temperature regime. Therefore, 
we conclude that the MSW method is useful as a com- 
plementary method to the ED analysis even in the pres- 
ence of moderate frustration as far as the ground state 
remains the ferromagnetic state. 

It is also predicted that the MSW approximation is 
reliable even at the limit of the stability of the ferro- 
magnetic ground state J\ = Ji s insofar as the exponent 
of the leading temperature dependence is concerned. Al- 
though we have explicitly demonstrated \ ~ T~ 4 ^ 3 on 
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havior should be observed even somewhat away from the 
transition points at finite temperatures. As the tempera- 
ture is lowered, the crossover to the true low temperature 
behavior in ferromagnetic or nonmagnetic phases should 
take place. Such a crossover behavior can be regarded as 
a precurser of the ground state phase transition if exper- 
imentally observed. 

In the present work, we have concentrated on the finite 
temperature properties in the parameter regime with fer- 
romagnetic ground state. However, our preliminary cal- 
culation suggests that the nonmagnetic phase consists of 
several different exotic phases with and without sponta- 
neous symmetry breakdown. The investigation of these 
phases will be reported elsewhere. 
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Appendix 

Let us define the density of states per site by 

w ^ = JL^ {6{x £a{k)) + 6{x £(,m} ■ (A ' 1} 

At low temperatures, only the gapless magnon mode a 
contributes. We consider the case where the dispersion 
relation of this gapless mode is given by e a (k) = Ak n . 
The density of states is then approximated as 



w 



2L ^ 2nnA- 



Rewriting (19) and (23) using the density of states, we 
have 



J o 



exp^T- 1 + v) - 1 2irnAi \n'J W ' 

(A-3) 



G/Mb) f°° exp(xT 1 + v)w(x)dx 



3T J (exp^rT- 1 + v) - 1) 



(5Mb) 2 T 



2 



3 2nnAn \nj \n 
where we define 

v = -p/T 



-I (-W- '-''l- (A4) 



(A-5) 



and F(a, v) is the Bose-Einstein integral function defined 

by 

F (a,v) = — ——■ (A-6) 



r(a) J e u+v - 1 

The behavior of this function for small v is known. 15 
For the present ourDOse. we onlv need the formula for 



nonintcger a. 

oo 

F(a, v) = T(l - a^"- 1 + J](/!)- 1 (-i;) i C(a - /)■ 

(A7) 



1=0 



For the calculation of the susceptibility, the partially in- 
tegrated form 

/ t — r T^du = T{a)F{a - 1, v) (A-8) 

J (e u+v - l) 2 v ' ' ' y ' 

is useful. The equations (A-3) and (A-4) can be expressed 
as 

(A-9) 



3 2-nnA^ \n 



*M a 4)' i " + <(H +0( ' 1: 

(A- 10) 

Solving (A-9) with respect to v and substituting into 
(A- 10), we find 

(<7Mb) Tl — 1 / _ . 1 . 7T\ 7T=T ^_ 

(2nSA*sm-) T 

3 n V nJ 

J 1 + iz^(L) T n\ n 



n - 1 \nj \nj 2nnSAn 



n(i-2n) c2 m r2 n 



2(n-l) 2 \nj \nj \2imSA 



1 



r Ti 



(All) 

Within the ferromagnetic phase, n = 2. Setting A = 
^ S, we have 

X 3T 2 ^ | 45^^^ ltorS*SS J " 

(A-12) 

At the stability limit of the ferromagnetic phase, n = 4. 
Hence, we have 

{ 9 ^?slAl2i \ 7 C(i)r(|)ri 

a — /i \ J- i 



Ti 



3 87rS*,43 



7/ c(|)r(i) . 1 1 

9 ^ SirSAi J j K ' 

In the low temperature limit, we have 

vT 4 / 3 Ai 

lim J = (A- 14) 

t^o(. 9Mb ) 2 2 V ; 

for 5=1/2. Numerically A ~ 0.04634074 for J r = 

xT 4/3 



J d = 1,(5 = 0.8, Ji = J u 



(3^b) 2 



0.1795937. This 



value is used in the fitting of the MS W0 data in Fig. 5. 
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